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' Abstract. In this paper, the problem of pinning control for synchronization of com- 

I plex dynamical networks is discussed. A cost function of the controlled network is 

\ defined by the feedback gain and the coupling strength of the network. An interest- 

, ing result is that lower cost is achieved by the control scheme of pinning nodes with 

ly^ , smaller degrees. Some rigorous mathematical analysis is presented for achieving 

lower cost in the synchronization of different star-shaped networks. Numerical sim- 
ulations on some non-regular complex networks generated by the Barabasi-Albert 
model and various star-shaped networks are shown for verification and illustration. 
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1 Introduction and problem formulation 

Complex networks are currently being studied across many fields of sciences, including physics, 
chemistry, biology, mathematics, sociology and engineering [H [21 [3l O [9l [151 [13 US]- A complex 
I network is a large set of interconnected nodes, in which a node is a fundamental unit with specific 
contents. Examples of complex networks include the Internet, food webs, cellular neural networks, 
^ \ biological neural networks, electrical power grids, telephone cell graphs, etc. Recently, synchro- 
nization of complex networks of dynamical systems has received a great deal of attention from the 
nonlinear dynamics community |101 112 1 [T6 l 1181 120j . A special control strategy called pinning control 
is used to achieve synchronization of complex networks; that is, only a fraction of the nodes or even 
a single node is controlled over the whole network [H [6l [TTl [2T] . This control method has become a 
common technique for control, stabilization and synchronization of coupled dynamical systems. In 
general, different nodes have different degrees in a network, thus a natural question is how different 
the effect would be when nodes with different degrees are pinned. 

Consider a dynamical network consisting of N identical and diffusively coupled nodes, with each 
node being an n-dimensional dynamical system. The state equations of the network are 

N 

Xi{t) = f{xi{t),t) + cJ2aijTxjit), i = l,2,---,N, (1) 
i=i 

where /(•) is the dynamical function of an isolated node, Xi = {xii,Xi2, ■ ■ ■ ,Xin) G -R" are the 
state variables of node i, constant c > represents the coupling strength, and T G ji^xN -g 
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inner linking matrix. Moreover, the coupling matrix A = {aij) E ^wx^v represents the coupling 

configuration of the network: If there is a connection between node i and node j {i ^ j), then 

aij = Qji = 1; otherwise, Uij = aji = (i ^ j); the diagonal entries of A are defined by 

JV 

aii = -'^aij, i = 1,2, • • • ,iV. (2) 

i=i 

Suppose that the network is connected in the sense of having no isolated clusters. Then, the 
coupling matrix A is irreducible. From Lemma 2 of |20J, it can be proved that zero is an eigenvalue 
of A with multiplicity one and all the other eigenvalues of A are strictly negative. 

Network ([1]) is said to achieve (asymptotical) synchronization if 

xi{t) X2{t) ^ • • • — > xj\f{t) s{t), as t — i- oo, (3) 

where, because of the diffusive coupling configuration, s{t) is a solution of an isolated node, which 
can be an equilibrium, a periodic or a chaotic orbit. As shown in [H El [HI EI] , this can be achieved 
by controlling several nodes (or even only one node) of the network. Without loss of generality, 
suppose that the controllers are added on the last N — k nodes of the network, so that the equations 
of the controlled network can be written as 

Xi{t) = f{xi{t),t) + cJ2f=i aijVxjit), z = 1,2, • • • ,k, 

(4) 

Xi{t) = f{xi{t),t) + cEf=i aijTxjit) - ceiT{xi{t) - s{t)), i = k + I, k + 2, ■ ■ ■ , N, 
where the feedback gains are positive constants. It can be seen that synchronizing all states Xi{t) 
to s{t) is determined by the dynamics of an isolated node, the coupling strength c > 0, the inner 
linking matrix F, the feedback gains > 0, and the coupling matrix A. 

As discussed in [6l lll[ [21\ . to achieve synchronization of complex dynamical networks, the con- 
trollers are generally preferred to be added to the nodes with larger degrees. However, it is also 
known that, to achieve a certain synchronizability of the network, the feedback gains £i usually 
have to be quite large. In [4j, when a single controller is used, the coupling strength c has to be 
quite large in general. From the view point of realistic applications, these are not expected and 
sometimes cannot be realized. Practically, a designed control strategy is expected to be effective 
and also easily implementable. In this paper, for various star-shaped networks and non-regular 
complex networks, a new concept of cost function is introduced to evaluate the efficiency of the 
designed controllers. It is found that surprisingly the cost can be much lower by controlling nodes 
with smaller degrees than controlling nodes with larger degrees. As will be seen, moreover, both the 
feedback gains Ei and the coupling strength c can be much smaller than those used in [U [6l [TTl [21] . 

The outline of this paper is as follows. In Section 2, a new definition of cost function and 
some mathematical preliminaries are given. Stability of different star-shaped networks controlled 
by pinning some nodes with small degrees are analyzed in Sections 3 and 4, respectively, where 
some simulated examples of dynamical networks are compared for illustration and verification. 
In Section 5, pinning control of non-regular complex dynamical networks of chaotic oscillators is 
studied through numerical simulations. Finally, Section 6 concludes the paper. 



2 The cost of pinning control 



Denote ej(t) = Xi{t) — s{t), where s{t) satisfies s{t) = f{s{t)). Then, the error equations of 
network ([T|) can be written as 

N 

e^{t) = f{xi{t),t) - f{s{t)) + cY.aij'^ejit), i = l,2,--- ,N, 

while the error equations of the controlled network (jj]) can be written as 

N 

e^{t) = f{xi{t),t)-f{s{t)) + c^dijTej{t), i = l,2,---,N, (5) 

i=i 

where dii = aii — ei,ei>0,i = k + l,k + 2,--- , N, and dij = aij otherwise. Let A = (dij) G R^^^ , 
and denote e{t) = (ei(t), e2{t), • • • , eAr(t))^. 
Differentiating ([5]) along s{t) gives 

e(t) = Difismeit) + cTeit)!^. (6) 

By analyzing the matrix A^, it is easy to see that all the eigenvalues of are negative, which are 
denoted by 

> Ai > • • • > Atv. 

There exists an orthogonal matrix U such that A^ = UJU"^, where J = diag{Xi, X2, ■ ■ ■ ,A7v}- 
Let e{t) = e{t)U. Then, from ([6]), one has 

^^{t) = [Df{sit)) + cX,r]ei{t), i = l,---,N. (7) 

Therefore, the local stability problem of network @ is converted into the stability problem of the 
independent linear systems ([7]). When the system function of an isolated node and the inner 
linking matrix T are fixed, the stability problem of systems ([7]) are dependent on the coupling 
strength c and the eigenvalues of A. Clearly, the smaller the eigenvalues of the matrix A are 
(Aj < 0, i = 1, • • • , A^), the smaller the coupling strength c > is needed to guarantee the same 
synchronizability of the network {[¥]), if systems d?]) have unbounded synchronization regions |18j . 
In the following, a cost function is introduced to describe the efficiency of the controllers. 

Definition 1 (Cost Function) Suppose that the feedback gain matrix is G = diag{ei, • • • , £n}, 
where ei>0,i = l,---,A^, are given as in The Cost Function is defined as 

N 

1=1 

Remark 1 The smaller the CF, the more efficient a control strategy to achieve the same goal 
of control, and the easier to be implemented. 

In order to discuss the effects of pinning control, the following Lemmas are needed. 

Lemma 2 [8j Let A = [aij] S (jnxn Hermitian, and let ann < • • • < ^22 ^ ^ii be a 
rearrangement of its diagonal entries in increasing order. Let the eigenvalues of A be ordered as 

Amin = A„ < Xfi-l < ■ • • < A2 < Ai = Amax- (8) 

Then 

(i) A„ < an < Ai for all i = 1, • • • , n, 

(ii) ail + 022 < A2, if Ai = 0. 



Remark 2 For diffusive networks, the smaller the A2, the easier the synchronization, if the 
network has an unbound synchronized region [18j. However, Lemma 2 shows that A2 is related to 
«ii + 022) where on and 022 are determined by two smallest nodes. Therefore, in this case, in order 
to improve the synchronizability, these small nodes should be pinned. 

Lemma 3 Let A G C'^x" be Hermitian, and 

i = [ ^1 A12 
~ [A^^ A2 + D \ ' 

where Ai G R'''"', Au G Rkx{N-k)^ ^ j^{N^k)x{N-k) jj ^ diag{ek+i, ■ ■ ■ ,eiv}- Then, 
A < —al, a > 0, if and only if Ai < —ali and A2 + D — Aj2{Ai + ali)~^ A12 < —C(l2, where Ii 
and I2 are identity matrixes with appropriate dimensions. 

Remark 3 Matrix D, i.e. £{ can be determined by the LMI method [7]. Suppose the eigenvalues 
of Ai are /3i > /32 > • • • > Pk- Then —a > (3i; and if £i —00, then —a j3i. 



3 Simple star-shaped networks 

Now, consider a network consisting of N identical nodes with a simple star-shaped coupling 
configuration: there exists a central node which connects to all the other non-central nodes, and 
there are no direct connections among the non-central nodes. Fig. 1 gives an example of such a 
network. 




Fig. 1 iV = 9 

The coupling matrix A of network ([T|) with a simple star-shaped coupling configuration can been 
written as 

" Ai A12 
A21 A2 
1) Gi?(^-i), ^21 =^r2, and A2 



A 



(9) 



'(Af-l)x(Af-l)- 



where Ai = -N + 1, A12 = (1, 1, ■ 

Let the feedback gain matrix be G = diag{Q,e, - ■ ■ ,e} G R^^^; that is, the controllers are 
added to all the non-central nodes for simplicity in this discussion. Then, the coupling matrix of 
the controlled network ^ is 

Ai A12 
A21 A2 - el 



A 



Theorem 1 If there are constants e > and k > 0, such that N — k — 1>0, e>k — 1 and 



e > 



k{N-k) 



N-k~l 

matrix A. 



then Xi{A) < —k, i = 1, - ■ ■ , N, where Xi{A), i = 1, - ■ ■ ,N, are the eigenvalues of the 



Proof: Consider the matrix 

A + kl 



Ai + k Ai2 
A21 {-l-e + k)I 

By the assumptions of the theorem, {—l — e + k)I<0. Let A = {Ai + kI) — Ai2{{—l — e + k)I)~'^A2i. 
Then 



A = -N+l+k+ 
< -N+l+k+ 
= 0. 



N-l 



1+e-k 



N~k~l 



fi(iV-l) 



Prom the Schur complement Lemma, A + kl < 0, which leads to the assertion of the theorem. □ 
Remark 4 Theorem 1 shows that under the control strategy of adding controllers to all the 
non-central nodes, as the constant e > increases, the eigenvalues of the matrix A will approach 
— (A^ — 1). Further, suppose that N > 2 and take k = 1. Then, based on Theorem 1, to guarantee 
that Xi{A) < —1, i = 1, ■ ■ ■ , N , the feedback gain e only needs to satisfy 

N-l 
' > N^.- 

In this case, if the coupling strength of network ([9]) is c, then the cost function is CF = ce{N — 1). 
From (jlOp and through direct calculations, it is easy to see that under the control strategy of 
adding controllers to the non-central nodes, CF > ^(j^f^^-^ suffices to guarantee that Xi{A) < 
-l,i = l,--- ,N. 

Remark 5 Consider the control strategy that a single controller is added to the central node. 
In this case, the feedback gain matrix is given by G = diag{s,0, • • • ,0} G fi^xN g^j^^ _ 
Let Xi{A) represent the largest eigenvalue of matrix A. Then, based on Lemma 2, Xi{A) > —1 
no matter how large the e, or equally the cost function CF, is taken. This shows that in order to 
guarantee Aj be small, the small nodes should be pinned. 

Some simulations on the network ([9]), shown in Fig. 1, are presented in Figs. 2~3, where the 
dynamics of an isolated node is a chaotic oscillator [T3], F = diag{0, 1, 0}, and the synchronized 
state is set to be s{t) = [7.9373 7.9373 21], which is an unstable equilibrium point of an isolated 
node 





(a) CF = c£ = 300c = 3000. (b) CF = ce{N - 1) = l.b{N - l)c = 120. 

Fig. 2 Synchronization of the simple star-shaped coupHng network ([9]) with couphng strength 
c — 10: (a) adding a single controller to the central node; (b) adding controllers to all the non-central nodes. 




(a) Ci^ = c£ = 500c = 3500. (b) CF = ce{N - 1) = l.b{N - l)c = 84. 



Fig. 3 Synchronization of the simple star-shaped coupling network ^ with coupling strength c = 7: 
(a) pinning the central node; (b) pinning all the non-central nodes. 

From Figs. 2~3, it can be seen that even with both a small coupling strength c and a small cost 
function CF, the synchronization of network ([9]) can be faster to achieve when the controllers are 
added to all the non-central nodes (nodes with smaller degrees) than the case that a single controller 
is added to the central node (the node with a larger degree). Although the more controllers are 
needed for pinning nodes with smaller degrees, the total cost is still lower and the control effect is 
better in comparison. 

4 Clusters of star-shaped networks with global coupling 

Consider a network consisting of identical nodes in clusters of star-shaped coupling configu- 
ration: there are k central nodes which are connected to each other; each central node may have 
different numbers of non-central nodes attached to it; there are no direct connections among the 
non-central nodes; and any central node has no direct connection to a non-central node attached 
to another central node. Fig. 4 gives an example of such a network for the case of A; = 3. 



Fig. 4 A cluster of star-shaped network with k = 3 
Without loss of generality, let the first k nodes be central nodes, each connecting to non- 
central nodes, where rii satisfy that < 77-2 < • • • < n^. Clearly, N = k + X]i=i ^i- '^^^ coupling 
matrix A of network ([1]) in such clusters of star-shaped coupling configuration can been written as 

Ai A12 
A21 A2 

where ^1 e ^12 S Rkx(N-k)^ = AJ2 G and A2 € RiN-k)x{N-k) _ ^. 

,1) € i = !,■■■ ,k. Then 



A 



A, 



-k + l-m 1 

1 -k + l-n2 



1 
1 

-A; + 1 - nfc 



12 



Ql 



q;2 



i21 



On 



Let the feedback gain matrix be G 
of the controlled network (jj]) is 







eL 



(N-k) 



A = A-G. 



. Then, the corresponding coupling matrix A 

(12) 



Theorem 2 If there are constants e > and A; > 0, such that ni > k, e > k—1 and e > M!ii+LJi), 
then Xi{A) < -k, i = 1, • • • ,N. 

Proof: Consider the matrix 

^ Ai + kl A12 

A21 {-l-e + k)I _ 

By the assumptions of the theorem, (-1-e+A;)/ < 0. Let ^ = (^i + A;/)-^i2((-l-e+A:)/)""^vl2i- 
Then A can be expressed as 



A + kl 



A 



Let nii 



-k + l-ni + k + -^ 

l+£—k 



1 

-k + l-n2 + k + 



n2 



l+e-k 



1 1 ■■■ -k + l-nk + k + - ^ 

-k + l-ni + k + + (A: - 1) • 1, i = l,2,--- ,k. Then 

rui = -ni + k + — ^ 
< -ni + ~k + ^rii 
= + k 

= ^(i-^)<o. 

According to the Gerschgorin disc theorem [8J, ^4 < 0. By the Schur complement Lemma, A + kl < 
0, which leads to the assertion of the theorem. □ 



Some simulations on the network (jlip shown in Fig. 4 are presented in Fig. 5, where the 
dynamics of an isolated node, the inner liking matrix F, and the synchronized state are the same 
as that given in Section 2. 




(a) CF = cek ^ SOO/cc = 9000. (b) CF = c£(X;Li "0 = 2.5(X;Li = 225. 

Fig. 5 Synchronization of network (jlip with coupling strength c — 10: (a) pinning k central nodes; 
(b) pinning all the non-central nodes. 

From Fig. 5, it can be seen that although all the non-central nodes are pinned, the cost is lower 
and the effect is better than the case where the central nodes are pinned. 

By using the Kronecker product, the controlled network ^ can be rewritten as 

e{t) = F{x,s)+c{A-G)®Te{t), (13) 

where e{t) = (ef(t), • • • , e%it)f, F(x, s) = - , • • • , (/(x^(t)) - fisit))ff, 

G = diag{0, ■ ■ ■ , 0, e, • • • , e} € i^^^^ and F = diag{ri,r2, • • • , r„}, = 1 or 0, i = 1, ■ • • , n. 

k N-k 

Theorem 3 If there are a positive diagonal matrix P and constants /Li > 0, e > 0, A; > 0, c > 0, 
such that rii > k, £ > k — 1, e > a,nd 

ni—k 

{x - y)^P{f{x, t) - f{y, t) - ckr{x - y)) < -^(x - y)^ {x - y), 

then the controlled network (jj]) is globally exponentially synchronized to s{t). 

Proof: Choose a Lyapunov function as V{t) = ^e^{t){I]\f (g) P)e{t). Its time derivative is 

Vit) = e^{t){lN^P)eit) 

= e^{t){lN (8) P){F{x, s) + c{A - G) (g) Te{t)) 

= e^{t){lN (8) P)F{x, s) + ce^{t){{A - G) (g, PT)e{t) 

= e^{t){lN (8) P){F{x, s) - c{kln ® F)e^(t)) + ce^ {t){{A - G) + kin) <8) Pre{t). 
By Theorem 2, {A - G) + kin < 0. So, 

and the theorem is thus proved. □ 



5 Some non-regular complex networks 



In this section, first, a non-regular coupled network consisting of 20 nodes is generated following 
the procedure of the well-known BA model. Fig. 6 shows the synchronization of the network with 
controllers being added to the three "biggest" nodes of degrees 15, 13 and 10, respectively. The 
dynamics of an isolated node, the inner liking matrix F, and the synchronized state are the same 
as that given in Section 2. From Figs. 6(c) and 6(d), it can be seen that the network does not 
synchronize any faster, although a larger feedback gain is used. In Fig. 7, a different control scheme 
is applied: eleven nodes with smaller degrees in the network are pinned with coupling strength c = 6 
and feedback gain e = 5, yielding cost function CF = 330. Compared with Fig. 6, a better control 
performance is obtained with a much smaller cost function. Figs. 8 and 9 also show that it is more 
efficient by pinning nodes with smaller degrees than pinning nodes with larger degrees. As shown 
in Fig. 8, to achieve a similar synchronization effect, a much smaller feedback gain is needed in 
pinning "smaller" nodes than pinning "larger" ones, thereby the cost function is also smaller in the 
former control scheme than that in the latter. Fig. 9 shows that although the coupling strength, 
the feedback gain and the cost function are similar, the synchronization effect is better in pinning 
"smaller" nodes than pinning "larger" ones. 




(c) (d) 



Fig. 6 Pinning the three "biggest" nodes with degrees 15, 13 and 10 in a 20-node non-regular 
coupled network: (a) c = 0, e 0, CF = 0. (b) c = 6, e = 0, CF = 0. (c) c = 6, e = 500, CF = 9000. (d) 
c = 6, e = 1000, CF = 18000. 
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Fig. 7 Pinning the eleven "smaller" nodes in a 20-node non-regular coupled network: c = 6, e = 5, 
CF = 330. 





(a) c = 8, e = 500, CF ^ 12000. 



(b) c = 6, £ = 8, CF = 528. 



Fig. 8 (a) Pinning the three "biggest" nodes with degrees 15, 13 and 10. (b) Pinning the eleven 
"smaller" nodes. 




(a) c = 6, £ = 55, CF = 660. (b) c = 6, e = 22, CF = 660. (c) c = 6, e = 10, CF = 660. 

Fig. 9 (a) Pinning the two "biggest" nodes with degrees 15 and 13. (b) Pinning the three "biggest" 
nodes with degrees 15, 13 and 10, and the two "sniaUest" nodes with the same degrees 3. (c) Pinning the 
eleven "smaller" nodes. 

It seems more efficient to use the control scheme of pinning nodes with smaller degrees than 
pinning nodes with larger degrees for achieving the desired synchronous states of some non-regular 
coupled dynamical networks. This is an interesting phenomenon that was not noticed before. 

6 Conclusions 

In this paper, pinning control to achieve synchronization of complex dynamical networks is 
further investigated. In contrary to the general perception, the nodes with smaller degrees play 
an important role in the synchronizability of some networks. It has been shown by computer 
simulations in this paper that to achieve a similar synchronization effect on the networks considered, 
a smaller coupling strength, a smaller feedback gain and thus a lower cost function are needed in 
the control scheme of pinning nodes with smaller degrees comparing with those needed in pinning 
nodes with larger degrees. In other words, it seems more efficient to use the control scheme 
of pinning nodes with smaller degrees than that of pinning nodes with larger degrees for some 
dynamical networks to achieve synchronization, an interesting phenomenon that deserves further 
investigation in the future. 
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1 Introduction and problem formulation 

Complex networks are currently being studied across many fields of sciences, including physics, 
chemistry, biology, mathematics, sociology and engineering [H [21 [3l O [9l [151 [13 US]- A complex 
I network is a large set of interconnected nodes, in which a node is a fundamental unit with specific 
contents. Examples of complex networks include the Internet, food webs, cellular neural networks, 
^ \ biological neural networks, electrical power grids, telephone cell graphs, etc. Recently, synchro- 
nization of complex networks of dynamical systems has received a great deal of attention from the 
nonlinear dynamics community |101 112 1 [T6 l 1181 120j . A special control strategy called pinning control 
is used to achieve synchronization of complex networks; that is, only a fraction of the nodes or even 
a single node is controlled over the whole network [H [6l [TTl [2T] . This control method has become a 
common technique for control, stabilization and synchronization of coupled dynamical systems. In 
general, different nodes have different degrees in a network, thus a natural question is how different 
the effect would be when nodes with different degrees are pinned. 

Consider a dynamical network consisting of N identical and diffusively coupled nodes, with each 
node being an n-dimensional dynamical system. The state equations of the network are 

N 

Xi{t) = f{xi{t),t) + cJ2aijTxjit), i = l,2,---,N, (1) 
i=i 

where /(•) is the dynamical function of an isolated node, Xi = {xii,Xi2, ■ ■ ■ ,Xin) G -R" are the 
state variables of node i, constant c > represents the coupling strength, and T G ji^xN -g 
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inner linking matrix. Moreover, the coupling matrix A = {aij) E ^wx^v represents the coupling 

configuration of the network: If there is a connection between node i and node j {i ^ j), then 

aij = Qji = 1; otherwise, Uij = aji = (i ^ j); the diagonal entries of A are defined by 

JV 

aii = -'^aij, i = 1,2, • • • ,iV. (2) 

i=i 

Suppose that the network is connected in the sense of having no isolated clusters. Then, the 
coupling matrix A is irreducible. From Lemma 2 of |20J, it can be proved that zero is an eigenvalue 
of A with multiplicity one and all the other eigenvalues of A are strictly negative. 

Network ([1]) is said to achieve (asymptotical) synchronization if 

xi{t) X2{t) ^ • • • — > xj\f{t) s{t), as t — i- oo, (3) 

where, because of the diffusive coupling configuration, s{t) is a solution of an isolated node, which 
can be an equilibrium, a periodic or a chaotic orbit. As shown in [H El [HI EI] , this can be achieved 
by controlling several nodes (or even only one node) of the network. Without loss of generality, 
suppose that the controllers are added on the last N — k nodes of the network, so that the equations 
of the controlled network can be written as 

Xi{t) = f{xi{t),t) + cJ2f=i aijVxjit), z = 1,2, • • • ,k, 

(4) 

Xi{t) = f{xi{t),t) + cEf=i aijTxjit) - ceiT{xi{t) - s{t)), i = k + I, k + 2, ■ ■ ■ , N, 
where the feedback gains are positive constants. It can be seen that synchronizing all states Xi{t) 
to s{t) is determined by the dynamics of an isolated node, the coupling strength c > 0, the inner 
linking matrix F, the feedback gains > 0, and the coupling matrix A. 

As discussed in [6l lll[ [21\ . to achieve synchronization of complex dynamical networks, the con- 
trollers are generally preferred to be added to the nodes with larger degrees. However, it is also 
known that, to achieve a certain synchronizability of the network, the feedback gains £i usually 
have to be quite large. In [4j, when a single controller is used, the coupling strength c has to be 
quite large in general. From the view point of realistic applications, these are not expected and 
sometimes cannot be realized. Practically, a designed control strategy is expected to be effective 
and also easily implementable. In this paper, for various star-shaped networks and non-regular 
complex networks, a new concept of cost function is introduced to evaluate the efficiency of the 
designed controllers. It is found that surprisingly the cost can be much lower by controlling nodes 
with smaller degrees than controlling nodes with larger degrees. As will be seen, moreover, both the 
feedback gains Ei and the coupling strength c can be much smaller than those used in [U [6l [TTl [21] . 

The outline of this paper is as follows. In Section 2, a new definition of cost function and some 
mathematical preliminaries are given. Stability of different star-shaped networks controlled by 
pinning some nodes with small degrees are analyzed in Sections 3, 4 and 5, respectively, where 
some simulated examples of dynamical networks are compared for illustration and verification. 
In Section 6, pinning control of non-regular complex dynamical networks of chaotic oscillators is 
studied through numerical simulations. Finally, Section 7 concludes the paper. 



2 The cost of pinning control 



Denote ej(t) = Xi{t) — s{t), where s{t) satisfies s{t) = f{s{t)). Then, the error equations of 
network ([T|) can be written as 

N 

e^{t) = f{xi{t),t) - f{s{t)) + cY.aij'^ejit), i = l,2,--- ,N, 

while the error equations of the controlled network (jj]) can be written as 

N 

e^{t) = f{xi{t),t)-f{s{t)) + c^dijTej{t), i = l,2,---,N, (5) 

i=i 

where dii = aii — ei,ei>0,i = k + l,k + 2,--- , N, and dij = aij otherwise. Let A = (dij) G R^^^ , 
and denote e{t) = (ei(t), e2{t), • • • , eAr(t))^. 
Differentiating ([5]) along s{t) gives 

e(t) = Difismeit) + cTeit)!^. (6) 

By analyzing the matrix A^, it is easy to see that all the eigenvalues of are negative, which are 
denoted by 

> Ai > • • • > Atv. 

There exists an orthogonal matrix U such that A^ = UJU"^, where J = diag{Xi, X2, ■ ■ ■ ,A7v}- 
Let e{t) = e{t)U. Then, from ([6]), one has 

^^{t) = [Df{sit)) + cX,r]ei{t), i = l,---,N. (7) 

Therefore, the local stability problem of network @ is converted into the stability problem of the 
independent linear systems ([7]). When the system function of an isolated node and the inner 
linking matrix T are fixed, the stability problem of systems ([7]) are dependent on the coupling 
strength c and the eigenvalues of A. Clearly, the smaller the eigenvalues of the matrix A are 
(Aj < 0, i = 1, • • • , A^), the smaller the coupling strength c > is needed to guarantee the same 
synchronizability of the network {[¥]), if systems d?]) have unbounded synchronization regions |18j . 
In the following, a cost function is introduced to describe the efficiency of the controllers. 

Definition 1 (Cost Function) Suppose that the feedback gain matrix is G = diag{ei, • • • , £n}, 
where ei>0,i = l,---,A^, are given as in The Cost Function is defined as 

N 

1=1 

Remark 1 The smaller the CF, the more efficient a control strategy to achieve the same goal 
of control, and the easier to be implemented. 

In order to discuss the effects of pinning control, the following Lemmas are needed. 

Lemma 2 [8j Let A = [aij] S (jnxn Hermitian, and let ann < • • • < ^22 ^ ^ii be a 
rearrangement of its diagonal entries in increasing order. Let the eigenvalues of A be ordered as 

Amin = A„ < Xfi-l < ■ • • < A2 < Ai = Amax- (8) 

Then 

(i) A„ < an < Ai for all i = 1, • • • , n, 

(ii) ail + 022 < A2, if Ai = 0. 



Remark 2 For diffusive networks, the smaller the A2, the easier the synchronization, if the 
network has an unbound synchronized region [18j. However, Lemma 2 shows that A2 is related to 
«ii + 022) where on and 022 are determined by two smallest nodes. Therefore, in this case, in order 
to improve the synchronizability, these small nodes should be pinned. 

Lemma 3 Let A G C'^x" be Hermitian, and 

i = [ ^1 A12 
~ [A^^ A2 + D \ ' 

where Ai G R'''"', Au G Rkx{N-k)^ ^ j^{N^k)x{N-k) jj ^ diag{ek+i, ■ ■ ■ ,eiv}- Then, 
A < —al, a > 0, if and only if Ai < —ali and A2 + D — Aj2{Ai + ali)~^ A12 < —C(l2, where Ii 
and I2 are identity matrixes with appropriate dimensions. 

Remark 3 Matrix D, i.e. £{ can be determined by the LMI method [7]. Suppose the eigenvalues 
of Ai are /3i > /32 > • • • > Pk- Then —a > (3i; and if £i —00, then —a j3i. 



3 Simple star-shaped networks 

Now, consider a network consisting of N identical nodes with a simple star-shaped coupling 
configuration: there exists a central node which connects to all the other non-central nodes, and 
there are no direct connections among the non-central nodes. Fig. 1 gives an example of such a 
network. 




Fig. 1 iV = 9 

The coupling matrix A of network ([T|) with a simple star-shaped coupling configuration can been 
written as 

" Ai A12 
A21 A2 
1) Gi?(^-i), ^21 =^r2, and A2 



A 



(9) 



'(Af-l)x(Af-l)- 



where Ai = -N + 1, A12 = (1, 1, ■ 

Let the feedback gain matrix be G = diag{Q,e, - ■ ■ ,e} G R^^^; that is, the controllers are 
added to all the non-central nodes for simplicity in this discussion. Then, the coupling matrix of 
the controlled network ^ is 

Ai A12 
A21 A2 - el 



A 



Theorem 1 If there are constants e > and k > 0, such that N — k — 1>0, e>k — 1 and 



e > 



k{N-k) 



N-k~l 

matrix A. 



then Xi{A) < —k, i = 1, - ■ ■ , N, where Xi{A), i = 1, - ■ ■ ,N, are the eigenvalues of the 



Proof: Consider the matrix 

A + kl 



Ai + k Ai2 
A21 {-l-e + k)I 

By the assumptions of the theorem, {—l — e + k)I<0. Let A = {Ai + kI) — Ai2{{—l — e + k)I)~'^A2i. 
Then 



A = -N+l+k+ 
< -N+l+k+ 
= 0. 



N-l 



1+e-k 



N~k~l 



fi(iV-l) 



Prom the Schur complement Lemma, A + kl < 0, which leads to the assertion of the theorem. □ 
Remark 4 Theorem 1 shows that under the control strategy of adding controllers to all the 
non-central nodes, as the constant e > increases, the eigenvalues of the matrix A will approach 
— (A^ — 1). Further, suppose that N > 2 and take k = 1. Then, based on Theorem 1, to guarantee 
that Xi{A) < —1, i = 1, ■ ■ ■ , N , the feedback gain e only needs to satisfy 

N-l 
' > N^.- 

In this case, if the coupling strength of network ([9]) is c, then the cost function is CF = ce{N — 1). 
From (jlOp and through direct calculations, it is easy to see that under the control strategy of 
adding controllers to the non-central nodes, CF > ^(j^f^^-^ suffices to guarantee that Xi{A) < 
-l,i = l,--- ,N. 

Remark 5 Consider the control strategy that a single controller is added to the central node. 
In this case, the feedback gain matrix is given by G = diag{s,0, • • • ,0} G fi^xN g^j^^ _ 
Let Xi{A) represent the largest eigenvalue of matrix A. Then, based on Lemma 2, Xi{A) > —1 
no matter how large the e, or equally the cost function CF, is taken. This shows that in order to 
guarantee Aj be small, the small nodes should be pinned. 

Some simulations on the network ([9]), shown in Fig. 1, are presented in Figs. 2~3, where the 
dynamics of an isolated node is a chaotic oscillator [T3], F = diag{0, 1, 0}, and the synchronized 
state is set to be s{t) = [7.9373 7.9373 21], which is an unstable equilibrium point of an isolated 
node 





(a) CF = c£ = 300c = 3000. (b) CF = ce{N - 1) = l.b{N - l)c = 120. 

Fig. 2 Synchronization of the simple star-shaped coupHng network ([9]) with couphng strength 
c — 10: (a) adding a single controller to the central node; (b) adding controllers to all the non-central nodes. 




(a) Ci^ = c£ = 500c = 3500. (b) CF = ce{N - 1) = l.b{N - l)c = 84. 



Fig. 3 Synchronization of the simple star-shaped coupling network ^ with coupling strength c = 7: 
(a) pinning the central node; (b) pinning all the non-central nodes. 

From Figs. 2~3, it can be seen that even with both a small coupling strength c and a small cost 
function CF, the synchronization of network ([9]) can be faster to achieve when the controllers are 
added to all the non-central nodes (nodes with smaller degrees) than the case that a single controller 
is added to the central node (the node with a larger degree). Although the more controllers are 
needed for pinning nodes with smaller degrees, the total cost is still lower and the control effect is 
better in comparison. 

4 Clusters of star-shaped networks with global coupling 

Consider a network consisting of identical nodes in clusters of star-shaped coupling configu- 
ration: there are k central nodes which are connected to each other; each central node may have 
different numbers of non-central nodes attached to it; there are no direct connections among the 
non-central nodes; and any central node has no direct connection to a non-central node attached 
to another central node. Fig. 4 gives an example of such a network for the case of A; = 3. 



Fig. 4 A cluster of star-shaped network with k = 3 
Without loss of generality, let the first k nodes be central nodes, each connecting to non- 
central nodes, where rii satisfy that < 77-2 < • • • < n^. Clearly, N = k + X]i=i ^i- '^^^ coupling 
matrix A of network ([1]) in such clusters of star-shaped coupling configuration can been written as 

Ai A12 
A21 A2 

where ^1 e ^12 S Rkx(N-k)^ = AJ2 G and A2 € RiN-k)x{N-k) _ ^. 

,1) € i = !,■■■ ,k. Then 



A 



A, 



-k + l-m 1 

1 -k + l-n2 



1 
1 

-A; + 1 - nfc 



12 



Ql 



q;2 



i21 



On 



Let the feedback gain matrix be G 
of the controlled network (jj]) is 







eL 



(N-k) 



A = A-G. 



. Then, the corresponding coupling matrix A 

(12) 



Theorem 2 If there are constants e > and A; > 0, such that ni > k, e > k—1 and e > M!ii+LJi), 
then Xi{A) < -k, i = 1, • • • ,N. 

Proof: Consider the matrix 

^ Ai + kl A12 

A21 {-l-e + k)I _ 

By the assumptions of the theorem, (-1-e+A;)/ < 0. Let ^ = (^i + A;/)-^i2((-l-e+A:)/)""^vl2i- 
Then A can be expressed as 



A + kl 



A 



Let nii 



-k + l-ni + k + -^ 

l+£—k 



1 

-k + l-n2 + k + 



n2 



l+e-k 



1 1 ■■■ -k + l-nk + k + - ^ 

-k + l-ni + k + + (A: - 1) • 1, i = l,2,--- ,k. Then 

rui = -ni + k + — ^ 
< -ni + ~k + ^rii 
= + k 

= ^(i-^)<o. 

According to the Gerschgorin disc theorem [8J, ^4 < 0. By the Schur complement Lemma, A + kl < 
0, which leads to the assertion of the theorem. □ 



Some simulations on the network (jlip shown in Fig. 4 are presented in Fig. 5, where the 
dynamics of an isolated node, the inner liking matrix F, and the synchronized state are the same 
as that given in Section 2. 




(a) CF = cek ^ SOO/cc = 9000. (b) CF = c£(X;Li "0 = 2.5(X;Li = 225. 

Fig. 5 Synchronization of network (jlip with coupling strength c — 10: (a) pinning k central nodes; 
(b) pinning all the non-central nodes. 

From Fig. 5, it can be seen that although all the non-central nodes are pinned, the cost is lower 
and the effect is better than the case where the central nodes are pinned. 

By using the Kronecker product, the controlled network ^ can be rewritten as 

e{t) = F{x,s)+c{A-G)®Te{t), (13) 

where e{t) = (ef(t), • • • , e%it)f, F(x, s) = - , • • • , (/(x^(t)) - fisit))ff, 

G = diag{0, ■ ■ ■ , 0, e, • • • , e} € i^^^^ and F = diag{ri,r2, • • • , r„}, = 1 or 0, i = 1, ■ • • , n. 

k N-k 

Theorem 3 If there are a positive diagonal matrix P and constants /Li > 0, e > 0, A; > 0, c > 0, 
such that rii > k, £ > k — 1, e > a,nd 

ni—k 

{x - y)^P{f{x, t) - f{y, t) - ckr{x - y)) < -^(x - y)^ {x - y), 

then the controlled network (jj]) is globally exponentially synchronized to s{t). 

Proof: Choose a Lyapunov function as V{t) = ^e^{t){I]\f (g) P)e{t). Its time derivative is 

Vit) = e^{t){lN^P)eit) 

= e^{t){lN (8) P){F{x, s) + c{A - G) (g) Te{t)) 

= e^{t){lN (8) P)F{x, s) + ce^{t){{A - G) (g, PT)e{t) 

= e^{t){lN (8) P){F{x, s) - c{kln ® F)e^(t)) + ce^ {t){{A - G) + kin) <8) Pre{t). 
By Theorem 2, {A - G) + kin < 0. So, 

and the theorem is thus proved. □ 



5 A network with multi-central nodes 



6 Some non-regular complex networks 

In this section, first, a non-regular coupled network consisting of 20 nodes is generated following 
the procedure of the well-known BA model. Fig. 8 shows the synchronization of the network 
with controllers being added to the three "biggest" nodes of degrees 15, 13 and 10, respectively. 
The dynamics of an isolated node, the inner liking matrix F, and the synchronized state are the 
same as that given in Section 2. From Figs. 8(c) and 8(d), it can be seen that the network does 
not synchronize any faster, although a larger feedback gain is used. In Fig. 9, a different control 
scheme is applied: eleven nodes with smaller degrees in the network are pinned with coupling 
strength c = 6 and feedback gain e = 5, yielding cost function CF = 330. Compared with Fig. 8, 
a better control performance is obtained with a much smaller cost function. Figs. 10 and 11 also 
show that it is more efficient by pinning nodes with smaller degrees than pinning nodes with larger 
degrees. As shown in Fig. 10, to achieve a similar synchronization effect, a much smaller feedback 
gain is needed in pinning "smaller" nodes than pinning "larger" ones, thereby the cost function is 
also smaller in the former control scheme than that in the latter. Fig. 11 shows that although the 
coupling strength, the feedback gain and the cost function are similar, the synchronization effect is 
better in pinning "smaller" nodes than pinning "larger" ones. 




(c) (d) 

Fig. 8 Pinning the three "biggest" nodes with degrees 15, 13 and 10 in a 20-node non-regular 
coupled network: (a) c = 0, e = 0, CF 0. (b) c = 6, e = 0, CF = 0. (c) c = 6, e = 500, CF = 9000. (d) 
c = 6, e = 1000, CF = 18000. 




0.5 1 



Fig. 9 Pinning the eleven "smaller" nodes in a 20-node non-regular coupled network: c = 6, e = 5, 
CF = 330. 




(a) c = 6, £ = 55, CF = 660. (b) c = 6, e = 22, CF = 660. (c) c = 6, e = 10, CF = 660. 

Fig. 11 (a) Pinning the two "biggest" nodes with degrees 15 and 13. (b) Pinning the three "biggest" 
nodes with degrees 15, 13 and 10, and the two "smaUest" nodes with the same degrees 3. (c) Pinning the 
eleven "smaher" nodes. 

It seems more efficient to use the control scheme of pinning nodes with smaller degrees than 
pinning nodes with larger degrees for achieving the desired synchronous states of some non-regular 
coupled dynamical networks. This is an interesting phenomenon that was not noticed before. 

7 Conclusions 

In this paper, pinning control to achieve synchronization of complex dynamical networks is 
further investigated. In contrary to the general perception, the nodes with smaller degrees play 
an important role in the synchronizability of some networks. It has been shown by computer 
simulations in this paper that to achieve a similar synchronization effect on the networks considered, 
a smaller coupling strength, a smaller feedback gain and thus a lower cost function are needed in 
the control scheme of pinning nodes with smaller degrees comparing with those needed in pinning 
nodes with larger degrees. In other words, it seems more efficient to use the control scheme 
of pinning nodes with smaller degrees than that of pinning nodes with larger degrees for some 
dynamical networks to achieve synchronization, an interesting phenomenon that deserves further 
investigation in the future. 
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